Inferring periodic orbits from spectra of simple shaped micro-lasers 
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Dielectric micro-cavities are widely used as laser resonators and characterizations of their spectra 
are of interest for various applications. We experimentally investigate micro-lasers of simple shapes 
(Fabry- Perot, square, pentagon, and disk). Their lasing spectra consist mainly of almost equidistant 
peaks and the distance between peaks reveals the length of a quantized periodic orbit. To measure 
this length with a good precision, it is necessary to take into account different sources of refractive 
index dispersion. Our experimental and numerical results agree with the superscar model describing 
the formation of long-lived states in polygonal cavities. The limitations of the two-dimensional 
approximation are briefly discussed in connection with micro-disks. 

PACS numbers: 42.55.Sa, 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

Two-dimensional micro-resonators and micro-lasers 
are being developed as building blocks for optical 
telecommunications [l|, [2[. Furthermore they are of in- 
terest as sensors for chemical or biological applications 
, [i| as well as billiard toy models for quantum chaos 
Towards fundamental and applied considerations, 
their spectrum is one of the main features. It was used, 
for instance, to experimentally recover some information 
about the refractive index [7[ or geometrical parameters 

a. 

In this paper we focus on cavities much larger than 
the wavelength and propose to account for spectra in 
terms of periodic orbit families. Cavities of the simplest 
and most currently used shapes were investigated: the 
Fabry-Perot resonator, polygonal cavities such as square 
and pentagon, and circular cavities. 

Our experiments are based on quasi two-dimensional 
organic micro-lasers [9|. The relatively straightfor- 
ward fabrication process ensures good quality and repro- 
ducibility as well as versatility in shapes and sizes (see 
Fig. [1]). The experimental and theoretical approaches 
developed in this paper can be easily extended to more 
complicated boundary shapes. Moreover this method is 
useful towards other kinds of micro-resonators, as it de- 
pends only on cavity shape and refractive index. 

The paper is organized as follows. In Section [III a de- 
scription of the two-dimensional model is provided to- 
gether with its advantages and limitations. In Section HU 
micro-lasers in the form of a long stripe are investigated 
as Fabry-Perot resonators to test the method and eval- 








FIG. 1: Optical microscope photographs of some organic 
micro-lasers: stripe (partial view, used as Fabry-Perot res- 
onator), square, pentagon, disk, quasi-stadium, and cardioid. 
Typical dimension: 100 \im. 



uate its experimental precision. This protocol is then 
further applied to polygonal cavities. In Section HV1 the 
case of square cavities is discussed whereas in Section |Vl 
dielectric pentagonal cavities are investigated. The the- 
oretical predictions based on a superscar model are com- 
pared to experiments as well as numerical simulations 
and a good agreement is found. Finally, in Section [Vj 
the case of several coexisting orbits is briefly dealt with 
on the example of circular cavities. 



II. PRELIMINARIES 
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Dielectric micro-cavities are quasi two-dimensional ob- 
jects whose thickness is of the order of the wavelength 
but with much bigger plane dimensions (see Fig. [1]). Al- 
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though such cavities have been investigated for a long 
time both with and without lasing, their theoretical de- 
scription is not quite satisfactory. In particular, the au- 
thors are not aware of true three-dimensional studies of 
high-excited electromagnetic fields even for passive cavi- 
ties. Usually one uses a two-dimensional approximation 
but its validity is not under control. 

Within such approximation fields inside the cavity and 
close to its two-dimensional boundary are treated differ- 
ently. In the bulk, one considers electromagnetic fields as 
propagating inside an infinite dielectric slab (gain layer) 
with refractive index n g \ surrounded by medias with re- 
fractive indices n\ and n 2 smaller than n g \. In our ex- 
periments, the gain layer is made of a polymer (PMMA) 
doped with a laser dye (DCM) and sandwiched between 
the air and a polymer (SOG) layer (see Fig. [2] (a) and [9[). 
It is well known (see e.g. [To| or [ll[) that in such ge- 
ometry there exist a finite number of propagating modes 
confined inside the slab by total internal reflection. The 
allowed values of transverse momentum inside the slab, 
g, are determined from the standard relation 
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where h is the slab thickness and T\^ 2 are the Fresnel re- 
flection coefficients on the two horizontal interfaces. For 
total internal reflection 
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FIG. 2: (a) Notations for refractive indexes and propagation 
wavenumbers. From top to bottom, the layers of our samples 
[9j] are the air (ni = 1), a polymer (PMMA) doped with a 
laser dye (DCM) (n g i — 1.54), and another polymer (SOG) 
(ri2 — 1.42) or silica (712 = 1.45). (b) Refractive index of the 
gain layer versus the wavelength inferred from ellipsometric 
measurements. 
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Here 6 is the angle between the direction of wave prop- 
agation inside the slab and the normal to the interface. 
The Vi parameter is 1 (resp. (n g i/rii) 2 ) when the mag- 
netic field (resp. the electric field) is perpendicular to 
the slab plane. The first and second cases correspond 
respectively to TE and TM polarizations. 

Denoting the longitudinal momentum, p = n g \k sin (9, 
as p — n e ffk, the effective refractive index, n e ff, is de- 
termined from the following dispersion relation 
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This equation has only a finite number of propagating 
solutions which can easily be obtained numerically. Fig. [3] 
presents possible propagating modes for our experimental 
setting m = 1 (air), n 2 = 1.42 (SOG) @ and n gX = 1.54 
deduced from ellipsometric measurements (see Fig.[2](b)) 
in the observation range. 



FIG. 3: Effective refractive index versus the thickness over 
wavelength variable calculated from Eq. (|4]). The refractive 
indices are assumed to be constant: 1 for air, 1.42 for SOG, 
and 1.54 for the gain layer (horizontal black lines). The TE 
polarization is plotted with solid blue lines and TM polariza- 
tion with dotted red lines. Integer / (see (HJ) increases from 
left to right starting from / = 0. 



The Maxwell equations for waves propagating inside 
the slab are thus reduced to the two-dimensional scalar 
Helmholtz equation: 



(A 



n 2 eff k 2 ) V in (x,y) = . 



(5) 



\I> represents the field perpendicular to the slab, i.e. the 
electric field for TM and the magnetic field for TE po- 
larization [29I ]. 

This equation adequately describes the wave propa- 
gation inside the cavity. But when one of these propa- 
gating modes hits the cavity boundary, it can partially 
escape from the cavity and partially be reflected inside 
it. To describe correctly different components of electro- 
magnetic fields near the boundary, the full solution of 
the three dimensional vectorial Maxwell equations is re- 
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quired, which to the authors knowledge has not yet been 
adressed in this context. Even the much simpler case of 
scalar scattering by a half-plane plate with a small but 
finite thickness is reduced only to numerical solution of 
the Wiener-Hopf type equation [121 ]. 

To avoid these complications, one usually considers 
that the fields can be separated into TE and TM po- 
larization and obey the scalar Helmholtz equations ([5|) 



(A 
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with rtin is the n e f f effective index inferred from Eq. (|4|) 
and n ou t the refractive index of the surrounding media, 
usually air so n out = 1. This system of two-dimensional 
equations is closed by imposing the following boundary 
conditions 
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Here r indicates the direction normal to the boundary 
and v depends on the polarization. When the electric 
(resp. magnetic) field is perpendicular to the cavity 
plane, called TM polarization (resp. TE polarization), 
v in, out = 1 (resp. v in, out = V^ou*)- Notice that these 
definitions of v are not the same for horizontal and ver- 
tical interfaces. 

We consider this standard two-dimensional approach 
keeping in mind that waves propagating close to the 
boundary (whispering gallery modes) may deviate 
significantly from two-dimensional predictions. In 
particular leakage through the third dimension could 
modify the life-time estimation of quasi-stationary states. 

Our polymer cavities are doped with a laser dye and 
uniformly pumped one by one from above [9], so that 
the pumping process induces no mode selection. The 
complete description of such lasing cavities requires the 
solution of the non-linear Maxwell-Bloch equations (see 
e.g. pjj Q EH and references therein). For clarity, 
we accept here a simplified point of view (see e.g. [1[ 
Sect. 24) according to which true lasing modes can be 
represented as a linear combination of the passive modes 
which may lase (i.e. for which gain exceeds losses) 



lasing 
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From physical considerations, it is natural to consider the 
\I/ m modes as the quasi-stationary states of the passive 
cavity. Though this choice leads to well known difficulties 
(see e.g. [1]) it is widely noticed and accepted at least 
for modes with small losses (cf. [3, El, EH ) • 

For each individual lasing mode, the C m coefficients 
could be determined only after the solution of the full 
Maxwell-Bloch equations. But due to the statistical na- 
ture of fluorescence the lasing effect starts randomly and 
independently during each pump pulse. So it is quite nat- 
ural to average over many pump pulses. Then the mean 
spectrum exhibits peaks at frequencies of all possible las- 
ing modes. The experimental data studied in this paper 



are recorded after integration over 30 pump pulses and 
agree with this simple statistical model. More refined 
verifications are in progress. 



III. FABRY-PEROT RESONATOR 

The Fabry-Perot configuration is useful for the calibra- 
tion control of further spectral experiments due to the 
non ambiguous single periodic orbit family which sus- 
tains the laser effect. 

A long stripe can be considered to a good approxima- 
tion as a Fabry-Perot resonator. In fact the pumping 
area is very small compared to the length (see Fig. [H 
(a)) and the material is slightly absorbing, so that reflec- 
tions at far extremities can be neglected. Moreover the 
pumping area is larger than the width of the stripe, thus 
the gain is uniformly distributed over the section. For a 
Fabry-Perot cavity, the emission is expected along both 
= and 6 = tt directions (see Fig. |U(a) for notations). 
Fig- SI (b) shows that this directional emission is observed 
experimentally which confirms the validity of our set-up. 
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FIG. 4: (a) Diagram summarizing the main features of the 
Fabry-Perot experiment, (b) Detected intensity versus angle 
for a Fabry-Perot experiment. 

The experimental spectrum averaged over 30 pump 
pulses is made up of almost regularly spaced peaks (see 
Fig. [5] (a)) which is typically expected for a Fabry- 
Perot resonator. In fact, due to coherent effects, the k 
wavenumbers of quasi-bound states of a passive Fabry- 
Perot cavity are determined from the quantization con- 
dition along the only periodic orbit of L = 2W length as 
for ©: 



rV L k 



eff(k) 
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where r is the Fresnel reflection coefficient and n e ff is 
the effective refractive index (j4]). The solutions of this 
equation are complex numbers: the imaginary part cor- 
responds to the width of the resonance and the real part 
(called k m afterwards) gives the position of a peak in the 
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spectrum and verify 



L k m n e ff(k m ) = 2tt m , m G N . 
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With Sk m = km+i—k m assumed to be small, the distance 
between adjacent peaks is constrained by 



Sk m [n e ff(k m ) 
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the full effective refractive index. It is a sum over 
two terms: one corresponding to the phase velocity, 
n e ff(km), and the other one to the group velocity, 

k m dri Q^ f (km)- If nfuii is considered as a constant over 
the observation range, which is true with a good accu- 
racy, 5k can be retrieved from the experimental spec- 
trum. For instance, the Fourier transform of the spec- 
trum (intensity versus k) is made up of regularly spaced 
peaks (Fig. [3(b) inset), with the first one (indicated with 
an arrow) centered at the optical length (L nf u u) and the 
others at its harmonics. So the geometrical length of the 
periodic orbit can be experimentally inferred from the 
knowledge of n f u u which is independently determined as 
described below. For the Fabry-Perot resonator, the ge- 
ometrical length is known to be 2W, thus allowing to 
check the experimental precision. The relative statistical 
errors on the W width is estimated to be less than 3 %. 
The error bars in Fig. [5] (b) are related to the first peak 
width of the Fourier transform and are less than 5 % of 
the optical length. 

The full effective refractive index, rif u ii, is indepen- 
dently inferred from ellipsometric measurement (Fig. [2] 
(b)) and standard effective index derivation described in 
the previous Section. Depending on the parameter h/X 
(thickness over wavelength), one or several modes are al- 
lowed to propagate. Our samples are designed such as 
only one TE and TM modes exist with n e ff effective 
refractive index according to Eq. 

In Fig. [3] the refractive index of the gain layer, n g i, 
is assumed to be constant: n g \ = 1.54 in the middle of 
the experimental window, A varying from 0.58 to 0.65 
/im. From Eq. (|4j) a n e ff = 1.50 is obtained in the ob- 
servation range with a h = 0.6 /im thickness, and corre- 
sponds to the phase velocity term. The group velocity 



term k n 



dn 



dk 



(k m ) is made up of two dispersion contri- 



butions: one from the effective index (about 4 %) and 
the other from the gain medium (about 7 %). The de- 
pendance of n g \ with the wavelength is determined with 
the GES 5 SOPRA ellipsometer from a regression with 
the Winelli II software (correlation coefficient: 0.9988) 
and plotted on Fig.[2](b). Taking into account all contri- 
butions (that means calculating the effective refractive 
index with a dispersed n g i), the nf u u full effective re- 
fractive index is evaluated to be 1.645 ± 0.008 in the 
observation range. So the group velocity term made up 
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FIG. 5: (a) Experimental spectrum of a Fabry-Perot res- 
onator with W — 150 \ivn. (b) Optical length versus Fabry- 
Perot width W. The experiments (red points) are linearly 
fitted by the solid red line. The dashed blue line corresponds 
to the theoretical prediction without any adjusted parameter. 
Inset: Normalized Fourier transform of the spectrum in (a) 
expressed as intensity versus wavenumber. 



of the two types of refractive index dispersion contribute 
for 10 % to the full effective index, which is significant 
compared to our experimental precision. The nj u u in- 
dex depends only smoothly on polarization (TE or TM), 
and on the h thickness, which is measured with a surface 
profllometer Veeco (Dektak 3 ST). Thus, the samples are 
designed with thickness 0.6 /im and the precision is re- 
ported on the full effective index which is assumed to be 
1.64 with a relative precision of about 1 % throughout 
this work. 

Considering all of these parameters, we obtain a satis- 
factory agreement between measured and calculated op- 
tical lengths, which further improves when taking into ac- 
count several Fabry-Perot cavities with different widths 
as shown on Fig. [5] (b). The excellent reproducibil- 
ity (time to time and sample to sample) is an addi- 
tional confirmation of accuracy and validity. With these 
Fabry-Perot resonators, we have demonstrated a spectral 
method to recover the geometrical length of a periodic 
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FIG. 6: (a) Experimental spectrum of a square-shaped micro- 
laser of 135 fim side width, (b) Optical length versus a square 
side width. The experiments (red points) are linearly fitted 
by the solid red line. The dashed blue line corresponds to the 
theoretical prediction (diamond periodic orbit) without any 
adjusted parameter. Top inset: Two representations of the 
diamond periodic orbit. Bottom inset: Normalized Fourier 
transform of the spectrum in (a) expressed as intensity versus 
wavenumber. 



orbit which can now be confidently applied to different 
shapes of micro-cavities. 



IV. SQUARE MICRO-CAVITY 

In the context of this paper square-shaped micro- 
cavities present a double advantage. Firstly, they are 
increasingly used in optical telecommunications [2|, Il7j ]. 
Secondly, the precision and validity of the parameters 
used above can be tested independently since there is 
only one totally confined periodic orbit family. In fact 
the refractive index is quite low (about 1.5), so the di- 
amond (see Fig. [6] (b), top inset) is the only short- 
period orbit without refraction loss (i.e. all reflection 



angles at the boundary are larger than the critical angle 
Xc — arcsin(l/n) ~ 42°.) 

In a square-shaped cavity light escapes mainly at the 
corners due to diffraction. Thus the quality design of cor- 
ners is critical for the directionality of emission but not 
for the spectrum. Indeed for reasonably well designed 
squared micro-cavity (see Fig. [I]), no displacement of the 
spectrum peaks is detectable by changing the 6 obser- 
vation angle. The spectra used in this paper are thus 
recorded in the direction of maximal intensity. 

Fig. [6] (a) presents a typical spectrum of a square- 
shaped micro-cavity. The peaks are narrower than in 
the Fabry-Perot resonator spectrum, indicating a better 
confinement, as well as regularly spaced, revealing a sin- 
gle periodic orbit. Data processing is performed exactly 
as presented in the previous Section: for each cavity the 
Fourier transform of the spectrum is calculated (Fig. [6] 
(b), bottom inset) and the position of its first peak is 
located at the optical length. Fig. [6] (b) summarizes the 
results for about twenty different micro-squares, namely: 
the optical length inferred from the Fourier transform 
versus the a square side width. These experimental 
results are fitted by the solid red line. The dotted blue 
line corresponds to an a priori slope given by nf u u 
(1.64) times the geometrical length of the diamond 
periodic orbit (L = 2\/2a). The excellent agreement 
confirms that the diamond periodic orbit family provides 
a dominant contribution to the quantization of dielectric 
square resonator. 

This result is far from obvious as square dielectric cav- 
ities are not integrable. At first glance the observed dom- 
inance of one short-period orbit can be understood from 
general considerations based on trace formulae which are 
a standard tool in semiclassical quantization of closed 
multi-dimensional systems (see e.g. [18|, [l9| and refer- 
ences therein). In general trace formulae express the 
density of states (and other quantities as well) as a sum 
over classical periodic orbits. For two-dimensional closed 
cavities 

d(k) = ^2S(k- k n ) « J2 c P e ikLp ~ ifIp + c.c. (13) 

n p 

where k is the wavenumber and k n are the eigenvalues 
of a closed cavity. The summation on the right part is 
performed over all periodic orbits labeled by p. L p is the 
length of the p periodic orbit, fi p is a certain phase accu- 
mulated from reflection on boundaries and caustics, and 
amplitude c p can be computed from classical mechanics. 
In general for integrable and pseudo-integrable systems 
(e.g. polygonal billiards) classical periodic orbits form 
continuous periodic orbit families and in two dimensions 

c P ~ (14) 
V ^p 

where A p is the geometrical area covered by a periodic or- 
bit family (see the example of circular cavities in Section 

EfJ). 
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Non-classical contributions from diffractive orbits and 
different types of creeping waves (in particular, lateral 
waves [HI) are individually smaller by a certain power 
of 1/k and are negligible in semiclassical limit k — > oo 
compared to periodic orbits. 

There exist no true bound states for open systems. 
One can only compute the spectrum of complex eigen- 
frequencies of quasi- stationary states. The real parts 
of such eigenvalues give the positions of resonances and 
their imaginary part measure the losses due to the leak- 
age from the cavity. 

For such systems it is quite natural to assume that the 
density of quasi-stationary states 

1 Im(fc n ) 



d(k) 



7T ^ (k- Re(k n )) 2 + lm(k n 



(15) 



can be written in a form similar to (jT3j) but the contri- 
bution of each periodic orbit has to be multiplied by the 
product of all reflection coefficients along this orbit (as it 
was done in a slightly different problem in fl9l]) 



d(fc)«£\ 



N p 

u 



.0) 



lfx " + c.c. . 



(16) 



Here N p is the number of reflections at the boundary and 

rp^ is the value of reflection coefficient corresponding to 
the j th reflection for the p periodic orbit. 

When the incident angle is larger than the critical 
angle the modulus of the reflection coefficient equals 1 
(see Eq. (|2|)), but if a periodic orbit hits a piece of 
boundary with angle smaller than the critical angle, then 
\r p \ < 1 thus reducing the contribution of this orbit. 
Therefore, the dominant contribution to the trace for- 
mula for open dielectric cavities is given by short-period 
orbits (c.p oc 1/ ^^/Lp~) which are confined by total internal 
reflection. For a square cavity with n = 1.5 the diamond 
orbit is the only confined short-period orbit which ex- 
plains our experimentally observation of its dominance. 

Nevertheless, this reasoning is incomplete because the 
summation of contributions of one periodic orbit and its 
repetitions in polygonal cavities does not produce a com- 
plex pole which is the characteristics of quasi-stationary 
states. 

In order to better understand the situation, we have 
performed numerical simulations for passive square cav- 
ities in a two-dimensional approximation with TM po- 
larization (see Section HJ and [20]). Due to symmetries, 
the quasi-stationary eigenstates can be classified accord- 
ing to different parities with respect to the square diago- 
nals. In Fig. (a), the imaginary parts of wavenumbers 
are plotted versus their real part for states antisymmet- 
ric according to the diagonals (that means obeying the 
Dirichlet boundary conditions along the diagonals) and 
called here ( ) states. 

These quasi-stationary states are clearly organized in 
families. This effect is more pronounced when wave func- 
tions corresponding to each family are calculated. For in- 
stance, wave functions for the three lowest families with 
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FIG. 7: (a) Imaginary parts versus real parts of the wavenum- 
bers of quasi- stationary states with (— — ) symmetry for a 
dielectric square resonator with n e // = 1.5 surrounded by air 
with n = 1. (b) The same as in (a) but for the states with the 
smallest modulus of the imaginary part (the most confined 
states). Inset. Empty triangles: the difference ([20]) between 
the real part of these wavenumbers and the asymptotic ex- 
pression. Filled circles: the same but when the correction 
term f[2T|) is taken into account. 



( ) symmetry are presented in Fig. [HI The other mem- 
bers of these families have similar patterns. The existence 
of such families was firstly noted in [21] for hexagonal di- 
electric cavities, then further detailed in [20[. 

One can argue that the origin of such families is anal- 
ogous to the formation of superscar states in pseudo- 
integrable billiards discussed in [22[ and observed ex- 
perimentally in microwave experiments in (23[. In gen- 
eral, periodic orbits of polygonal cavities form continuous 
families which can be considered as propagating inside 
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FIG. 8: Squared modulus of wave functions with 
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symmetry calculated with numerical simulations, (a) ak = 68.74— .026 : 
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(a) (b) (c) 

FIG. 9: Squared modulus of wave functions calculated within the superscar model ([17} and corresponding to the parameters 
of Fig. E] 




FIG. 10: Unfolding of the diamond periodic orbit. Thick lines 
indicate the initial triangle. 



straight channels obtained by unfolding classical motion 
(see Fig. [TO]) . These channels (hatched area Fig. [TQ|) are 
restricted by straight lines passing through cavity cor- 
ners. In [22] it was demonstrated that strong quantum 
mechanical diffraction on these singular corners forces 
wave functions in the semiclassical limit to obey simple 



boundary conditions on these (fictitious) channel bound- 
aries. More precisely it was shown that for billiard prob- 
lems \I/ on these boundaries take values of the order of 
0(1/Vk) — > when k — > oo. This result was obtained 
by using the exact solution for the scattering on peri- 
odic array of half-planes. No such results are known for 
dielectric problems. Nevertheless, it seems natural from 
semiclassical considerations that a similar phenomenon 
should appear for dielectric polygonal cavities as well. 

Within such framework, a superscar state can be con- 
structed explicitly as follows. After unfolding (see Fig. 
[TO]) , a periodic orbit channel has the form of a rectan- 
gle. Its length equals the periodic orbit length and its 
width is determined by the positions of the closest sin- 
gular corners. The unfolded superscar state corresponds 
to a simple plane wave propagating inside the rectangle 
taking into account all phase changes. It cancels at the 
fictitious boundaries parallel to the x direction and is 
periodic along this direction with a periodicity imposed 
by the chosen symmetry class. This procedure sets the 
wavenumber of the state and the true wavefunction is 
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obtained by folding back this superscar state. 

Superscar wave functions with ( ) symmetry asso- 
ciated with the diamond orbit (see Fig. [10]) are expressed 
as follows: 



+ sin (ftm ^x — 28^ sin ^ 



jW') (17) 

where x' and y' are coordinates symmetric with respect 
to square side. In coordinates as in Fig. [10] 

x' = y , ?/ = z 

In ([TT]) m and p are integers with p = 1,2,..., and 
m ^> 1. Z = \/2a is the half of the diamond periodic 
orbit length [30[, J is the phase of the reflection coeffi- 
cient defined by r = exp(— 2iS). For simplicity, we ig- 
nore slight changes of the reflection coefficient for differ- 
ent plane waves in the functions above. So S is given by 
([3]) with v = 1 for TM polarization and 9 = 7r/4. And 

ttm ^ is the momentum defined by 



2nm 



(18) 



This construction conducts to the following expression 
for the real part of the wavenumbers [31] 



27r^/(ra+ -S) 2 



- p* 



27r(ra+ -5) + 0( — ) . (19) 
7r m 



To check the accuracy of the above formulae we plot in 
Fig. [9] scar wave functions ([TT]) with the same parameters 
as those in Fig. [8] The latter were computed numerically 
by direct solving the Helmholtz equations ([6]) but the 
former looks very similar which supports the validity of 
the superscar model. 

The real part of the wavenumbers is tested too. In 
Fig. [7] (b) the lowest loss states (with the smallest mod- 
ulus of the imaginary part) with (— — ) symmetry are 
presented over a larger interval than in Fig. [T] (a). The 
real parts of these states are compared to superscar pre- 
dictions ([T9]) with p = 1, leading to a good agreement. 
To detect small deviations from the theoretical formula, 
we plot in the inset of Fig. [T] (b) the difference between 
a quantity inferred from numerical simulations and its 
superscar prediction from ([T9]) . 



SE = ( ^Rek 



(m+-5) 2 +p< 

7T 



(20) 



From this curve it follows that this difference tends to 
zero with m increasing, thus confirming the existence of 
the term proportional to p 2 . By fitting this difference 
with the simplest expression 





(a) 



(b) 



SE 



c 

rn 



(21) 



FIG. 11: Simplest whispering gallery periodic orbit family 
for a pentagonal cavity, (a) Solid line indicates the inscribed 
pentagon which is an isolated periodic orbit. A periodic orbit 
in its vicinity is plotted with dashed line. It belongs to the 
family of the five-pointed star periodic orbit. The fundamen- 
tal domain is indicated in grey, (b) Boundary of the family 
of the five-pointed star periodic orbit. 



we find that c ~ —6.9. By subtracting this correction 
term from the difference ([20]) . one gets the curve indi- 
cated with filled circles in inset of Fig. [T] (b). The result 
is one order of magnitude smaller than the difference it- 
self. 

All these calculations confirm that the real parts of 
resonance wavenumbers for square dielectric cavities are 
well reproduced in the semiclassical limit by the above 
superscar formula ([T9]) and our experimental results can 
be considered as an implicit experimental confirmation 
of this statement. 



V. PENTAGONAL MICRO-CAVITY 

The trace formula and superscar model arguments can 
be generalized to all polygonal cavities. The pentagonal 
resonator provides a new interesting test. In fact, due to 
the odd number of sides, the inscribed pentagonal orbit 
(indicated by solid line in Fig. [TT] (a)) is isolated. The 
shortest confined periodic orbit family is twice longer. 
It is represented with a dashed line in Fig. [TT] (a) and 
can be mapped onto the five-pointed star orbit drawn in 
Fig. [TT] (b) by continuous deformation. In this Section 
we compare the predictions of the superscar model for 
this periodic orbit family with numerical simulations and 
experiments. 

Due to the C$ v symmetry, pentagonal cavities sustain 
10 symmetry classes corresponding to the rotations by 
2tt/5 and the inversion with respect to one of the sym- 
metry axis. We have studied numerically one symmetry 
class in which wave functions obey the Dirichlet bound- 
ary conditions along two sides of a right triangle with 
angle 7r/5 (see Fig. [TT] (a) in grey). The results of these 
computations are presented in Fig. [12] 

As for the square cavity, lowest loss states are orga- 
nized in families. The wave functions of the three lowest 
loss families are plotted in Fig. [13] and their superscar 
structure is obvious. 

The computation of pure superscar states can be per- 
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outside the periodic orbit channel. 

The quantized values of the longitudinal momentum, 
ft, are obtained by imposing that the function ([22]) is 
periodic along the channel when all phases due to the 
reflection with the cavity boundaries are taken into ac- 
count 



kL = 2?r ( M 



10 



(23) 



Here M is an integer and L is the total periodic orbit 
length. For the five-pointed star orbit (see Fig. [TT]) 



L — 10a cos( — ) 
5 



(24) 



and S is the phase of the reflection coefficient given by 
® with v = 1 (for TM polarization) and = 3tt/10. 
For these states the real part of the wavenumber is the 
following 



nLRek = 2tt M 



Wave function inside the cavity are obtained by folding 
back the scar function ([22]) and choosing the correct rep- 
resentative of the chosen symmetry class. When Dirich- 
let boundary conditions are imposed along two sides of a 
right triangle passing through the center of the pentagon 
(see Fig. [15]), M must be written as M = 5 (2m) if p is 
odd and M = 5 (2m — 1) if p is even. Then the wave 
function inside the triangle is the sum of two terms 



(b) 



FIG. 12: (a) Wave numbers for a pentagonal cavity, a is the 
side length of the cavity. The three most confined families 
are indicated by solid, dashed and long-dashed lines, (b) The 
difference (|28|) between the real part of quasi-energies and 
superscar expression ([25]) for the three indicated families in 
(a). 



formed as in the previous Section. The five-pointed star 
periodic orbit channel is shown in Fig. [151 In this case 
boundary conditions along horizontal boundaries of pe- 
riodic orbit channel are not known. By analogy with su- 
perscar formation in polygonal billiards [22[, we impose 
that wave functions tend to zero along these boundaries 
when k — > oo. 

Therefore, a superscar wave function propagating in- 
side this channel takes the form 

7T 

^scar(^, y) = exp(ira) sin(— py)Q(y)Q(w - y) , (22) 

w 

where w is the width of the channel (for the five-pointed 
star orbit w = a sin(7r/5) where a is the length of the pen- 
tagon side). Q(x) is the Heavyside function introduced 
here to stress that superscar functions are zero (or small) 



^m iP (x, y) = sin(ft m x) sin(— py)Q(y)Q(w - y) + 



w 



+ sin(^ m x / — 25) sm(—py')Q(y')Q(w — y') 



where the longitudinal momentum is 
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27r(m+ -&-£) 



(26) 



(27) 



with £ = for odd p and £ = 1/2 for even p. x' and y' 
in ([26]) are coordinates of the point symmetric of (x, y) 
with respect to the inversion on the edge of the pentagon. 
In the coordinate system when the pentagon edge passes 
through the origin (as in Fig. [T5j) 

x r = x cos 2<j)-\-y sin 20 , y' = x sin 20 — y cos 20 

and = 7r/5 is the inclination angle of the pentagon 
side with respect to the abscissa axis. Wavefunctions ob- 
tained with this construction are presented in Fig. [Ml 
They correspond to the first, second, and third perpen- 
dicular excitations of the five-star periodic orbit family 
{p = 1, 2, and 3). 

To check the agreement between numerically computed 
real parts of the wavenumbers and the superscar predic- 
tion ([25]) and (fT7|). we plot in Fig. [12] (b) the following 
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(a) (b) (c) 

FIG. 14: Squared modulus of wave functions calculated within the superscar model and corresponding to the parameters of 
Fig. m 




FIG. 15: Unfolding of the five- star periodic orbit for a pen- 
tagonal cavity. Thick lines indicate the initial triangle. 



scales in Figs. 171 and [T2]) , we found it convenient to fit 
numerically the £ constant. The best fit gives £ ~ 0.44, 
2.33, and 5.51 for the three most confined families (for 
pure scar functions this constant is 1, 4, 9 respectively). 
The agreement is quite good with a relative accuracy 
of the order of 10 -4 (see Fig. [12] (b)). Irrespective of 
precise value of £ the total optical length, nL, is given 
by ([25]) and leads to an experimental prediction twice 
longer than the optical length of the inscribed pentagon, 
which is an isolated periodic orbit and thus can not base 
superscar wavefunctions. 



difference 
5E 



ncos( 
2tt 



-Keka ri 



rn H — S — £ 

7T 



2 tan( 



• (28) 



" 2 As our numerical sim- 



For pure scar states £ = p" . 
ulations have not reached the semiclassical limit (see 



Comparison with experiments confirms the superscar 
nature of the most confined states for pentagonal res- 
onators. In fact, the spectrum and its Fourier transform 
in Fig. [16] correspond to a pentagonal micro-laser with 
side a = 80 /im, and show a periodic orbit with optical 
length 1040 ± 30 jim to be compared with the five-star 
optical length nf u ul0a cos(7r/5) = 1061 jam. The agree- 
ment is better than 2%. 

This result is reproducible for cavities with the same size. 
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Other sizes have been tested as well. For smaller cavities, 
the five-pointed star orbit is not identifiable due to lack 
of gain, whereas for bigger ones it is visible but mixed 
with non confined periodic orbits. This effect, not spe- 
cific to pentagons, can be assigned to the contribution of 
different periodic orbit families which become important 
when the lasing gain exceeds the refractive losses. We will 
describe this phenomenon in a future publication [23 ]. 

The good agreement of numerical simulations and ex- 
periments with superscar predictions gives an additional 
credit to the validity of this approach even for non-trivial 
configurations. 




FIG. 17: (a) Two examples of periodic orbits: the square and 
the pentagon, (b) Two representations of the square periodic 
orbit and the caustic of this family in red. 
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FIG. 16: Experimental spectrum of a pentagonal micro-laser 
of 80 fim side length. Inset: Normalized Fourier transform of 
the spectrum expressed as intensity versus wavenumber. 



VI. MICRO-DISKS 

Micro-disk cavities are the simplest and most widely 
used micro-resonators. In the context of this work, they 
are of interest because of the coexistence of several pe- 
riodic orbit families with close lengths. For low index 
cavities (n ~ 1.5) each regular polygon trajectory with 
more than four sides is confined by total internal reflec- 
tion. 

In the two-dimensional approximation passive circu- 
lar cavities are integrable and the spectrum of quasi- 
stationary states can be computed from an explicit quan- 
tization condition 



l Jm(nkR) V H { rk\kR) 



(29) 



Here R is the radius of the disk, n the refractive index of 
the cavity, and v = 1 (resp. v = n 2 ) for the TM (resp. 
TE) polarization. For each angular quantum number m, 
an infinite sequence of solutions, /c m , g , is deduced from 
([29]) . They are labeled by the q radial quantum number. 



For large \k\ the fc m ,z wavenumbers are obtained from a 
semiclassical expression (see e.g. |25|) and the density of 
quasi-stationary states ([T5]) can be proved to be rewrit- 
ten as a sum over periodic orbit families. The derivation 
of this trace formula assumes only the semi-classical ap- 
proximation (\k\R^> 1) and can be done in a way similar 
to that of the billiard case (see e.g. [27[), leading to an 
expression closed to (fT6j) 



A 



d(k) oc —^= \r p \ Np cos(nL p k- 



(30) 



the p index specifies a periodic orbit family, 
formula depends on periodic orbit parame- 
the number of bounces on the boundary, N p , 



Here 
This 

ters: 

the incident angle on the boundary, Xp? the length, 
L p = 2N p Rcos(x P ), and the area covered by periodic 
orbit family, A p = ttR 2 cos 2 (xp), which is the area 
included between the caustic and the boundary (see Fig. 
im (b)). 2S P is the phase of the reflection coefficient at 
each bounce on the boundary (see Eq. ([3])) and \r p \ is 
its modulus. 

For orbits confined by total internal reflection S p does 
not depend on kR in the semi-classical limit, and r p is 
exponentially close to 1 [25|, [26[. From (|3Q|) it follows 
that each periodic orbit is singled out by a weighing 
coefficient c p = -A^= \r p \ Np . Considering the experi- 

mental values \k\R ~ 1000, \r p \ can be approximated 
to unity with a good accuracy for confined periodic 
orbits, and thus c n = -4= depends only on geometrical 



quantities. Fig. [TBI shows the evolution of c p for polygons 
when the number of sides is increasing. As the critical 
angle is close to 45°, the diameter and triangle periodic 
orbits are not confined and the dominating contribution 
comes from the square periodic orbit. So we can 
reasonably conclude that the spectrum ([T5]) of a passive 
two-dimensional micro-disk is dominated by the square 
periodic orbit. 

The experimental method described in the previous 
Sections has been applied to disk-shaped micro-cavities. 
A typical experimental spectrum is shown on Fig. [19] (a). 
The first peak of its Fourier transform (see Fig. [19] (b) 
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FIG. 18: Vertical red sticks: c p coefficient for polygons con- 
fined by total internal reflection (square, pentagon, hexagon, 
etc.). The dotted blue line indicates the position of the 
perimeter. 



inset) has a finite width coming from the experimental 
conditions (discretization, noise, etc..) and the contribu- 
tions of several periodic orbits. This width is represented 
as error bars on graph [19] (b). The continuous red line 
fitting the experimental data is surrounded by the dashed 
green line and the dotted blue line corresponding to the 
optical length of the square and hexagon respectively, 
calculated with nf u u = 1.64 as in the previous Sections. 
The perimeter (continuous black line) overlaps with a 
large part of the error bars which evidences its contribu- 
tion to the spectrum, but it is not close to experimental 
data. 

These experimental results seem in good agreement 
with the above theoretical predictions. But actu- 
ally these resonances, usually called whispering gallery 
modes, are living close to the boundary. Thus both 
roughness and three-dimensional effects must be taken 
into account. At this stage it is difficult to evaluate and 
to measure correctly such contributions for each periodic 
orbit. For micro-disks with a small thickness (about 0.4 
jam) and designed with lower roughness, the results are 
more or less similar to those presented in Fig. [19] (b). 
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FIG. 19: (a) Experimental spectrum of a micro-disk of 30 \ivn 
radius, (b) Optical length versus radius. The experiments 
(red points) are linearly fitted by the solid red line. The 
other lines correspond to theoretical predictions without any 
adjusted parameters: the dashed green line to the square, the 
dotted blue line to the hexagon, and the solid black line to 
the perimeter. Inset: Normalized Fourier transform of the 
spectrum in (a) expressed as intensity versus wavenumber. 



VII. CONCLUSION 

We demonstrate experimentally that the length of the 
dominant periodic orbit can be recovered from the spec- 
tra of micro-lasers with simple shapes. Taking into ac- 
count different dispersion corrections to the effective re- 
fractive index, a good agreement with theoretical predic- 
tions has been evidenced first for the Fabry-Perot res- 
onator. Then we have tested polygonal cavities both 
with experiments and numerical simulations, and a good 
agreement for the real parts of wavenumbers has been 
obtained even for the non trivial configuration of the pen- 
tagonal cavity. 

The observed dominance of confined short-period orbits 
is, in general, a consequence of the trace formula and the 



formation of long-lived states in polygonal cavities is re- 
lated to strong diffraction on cavity corners. 
Finally, the study of micro-disks highlights the case of 
several orbits and the influence of roughness and three- 
dimensional effect. 

Our study opens the way to a systematic exploration of 
spectral properties by varying the shape of the boundary. 
In increasing the experimental precision even tiny details 
of trace formulae will be accessible. The improvement of 
the etching quality will suppress the leakage due to sur- 
face roughness and lead to a measure of the diffractive 
mode losses which should depend on symmetry classes. 
From the point of view of technology, it will allow a better 
prediction of the resonator design depending on the ap- 
plications. From a more fundamental physics viewpoint, 
it may contribute to a better understanding of open di- 
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electric billiards. 
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